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Abstract. This paper deals with x^^\ the six-particle contribution to the 
magnetic susceptibility of the square lattice Ising model. We have generated, 
modulo a prime, series coefficients for x'^' ■ The length of the series is sufficient to 
produce the corresponding Fuchsian linear differential equation (modulo a prime). 
We obtain the Fuchsian linear differential equation that annihilates the "depleted" 
series il>(^) = ^(S) — ^X^*^ + ^x''^'- The factorization of the corresponding 
differential operator is performed using a method of factorization modulo a prime 
introduced in a previous paper. The "depleted" differential operator is shown 
to have a structure similar to the corresponding operator for x'^"*- H splits into 
factors of smaller orders, with the left-most factor of order six being equivalent to 
the symmetric fifth power of the linear differential operator corresponding to the 
elliptic integral E. The right-most factor has a direct sum structure, and using 
series calculated modulo several primes, all the factors in the direct sum have 
been reconstructed in exact arithmetics. 
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1. Introduction and recalls 

The magnetic susceptibility (high temperature x+ ^-nd low temperature X-) of the 
square lattice Ising model is given by [1] 
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and 

in terms of the self-dual temperature variable w — ^s/ (1 + s^), with s = sinh(2 J/fcT). 
The n-particle contributions x*-"-* are given by n — 1 dimensional integrals [11[31I11[S], 

'C")^ ;^-(n/ *)(n .,)■«•■■' -(<=-). (3) 



=1 



wher^l 



and 



with 



TJ h h - 2 sin (((/), -(/)j-)/2)-yx-x 



1 < i < J < n 



(4) 



^(„) ^ l±njU^^ (5) 



2u; 



1 - 2wcos(0,) + 2wcos{<p^)f - 4u;2 



(6) 



2u; 



^ ' ^-.. = (7) 

(1 - 2wcos(0i)) - 4u)2 j=i 

As n grows the series generation in the variable w of the integrals ^ becomes 
very time consuming. In [B] calculations modulo a prime were performed on simplified 

(n) 

integrals and this work demonstrated that most of the pertinent information 
(singularities, critical exponents, ...) can be obtained from linear ODEs known modulo 
a prime corresponding to the integrals ^^^^ . In order to go beyond x*-*^ this strategy 
was adopted previously for the 5-particle contribution x^^-' [71 [8] and here for the 
6-particle contribution x^' • 

In a previous paper massive computer calculations were performed on x*-^-* , x''^' 
and X (in exact arithmetics and/or modulo a prime). These calculations confirmed 
previously conjectured singularities for the linear ODEs of the x^"'''s as well as their 
critical exponents, and shed some light on important physical problems such as the 
existence of a natural boundary for the susceptibility of the square Ising model and 
the subtle resummation of logarithmic behaviours of the n-particle contributions x^""* 
to give rise to the power laws of the full susceptibility x- As far as x'^^'' is concerned, 
the linear ODE for x^*^^ was found modulo a single prime [7] and it is of minimal order 
33. 



1.1. Results on x^^^ 

In [8] the linear differential operator for x^^-* was carefully analysed. In particular it 
was found that the minimal order linear differential operator for x^^"* can be reduced to 
a minimal order linear differential operator L29 of order 29 for the linear combination 

§ The Fermionic term G^"^ has several representations [3]. 
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We shall use the term "depleted" series for a series obtained by substracting from 
X^"^ a definite amount of the lower n-particle contributions x^""^''-', k = 1, 2, • • •, as 
in ((SI), such that the differential operator annihilating the depleted series is of lower 
order. Since the depleted series is annihilated by an ODE of lower order, it follows 
that in the ODE for the original series, we must have the occurrence of a direct sum 
structure. It was found [8] that the linear differential operator L29, can be factorised 
as a product of an order five, an order twelve, an order one, and an order eleven linear 
differential operator 

L2g = L5 ■ L12 ■ Li ■ Lii, (9) 
where the order eleven linear differential operator has a direct-sum decomposition 

Ln = {Z2 ■ Ni) ® F2 © (F3 • F2 ■ LI) . (10) 
Z2 is a second order operator also occurring in the factorization of the linear 
differential operator [9] associated with x^^'^ and it corresponds to a modular form 
of weight one [10]. V2 is a second order operator equivalent to the second order 
operator associated with x*^^) (or equivalently to the complete elliptic integral E). 
F2 and F3 are remarkable second and third order globally nilpotent linear differential 
operators [H llOj . The first order linear differential operator Li quite remarkably 
has a polynomial solution. The fifth order linear differential operator L5 was shown 
to be equivalent to the symmetric fourth power of (the second order operator) Le 
corresponding to the complete elliptic integral E. The complete and detailed analysis 
of L12, the order twelve operator in © is beyond our current computional ressources 
(see [8] for details). 

It is important to note that these factorization results are exact and have been 
obtained from series and ODEs obtained modulo a single prime. For the reconstruction 
in exact arithmetics of the factors occurring in the differential operator in, we had 
to obtain the series and ODEs for more than one prime. The length of the series 
necessary to obtain the underlying ODE is initially unknown, except perhaps for some 
rough estimates. Once the first non-minimal order ODEs have been obtained modulo 
a prime, the minimum length of the series necessary to obtain non-minimal order 
ODEs for any other primes is known exactly. This knowledge comes from a relation 
we reported in [7] and that we called the " ODE formula" . Beyond understanding 
the terms occurring in the "ODE formula" and the light they shed on the ODEs 
underlying the problem, the formula has been of most importance in terms of gains in 
the computational effort. For instance, we initially generated, modulo a prime, 10000 
terms for x*^^"* and we found that we can obtain non-minimal ODEs using only some 
7400 terms, while non-minimal order ODEs for can be obtained using some 6200 
terms, representing a great reduction in the required computational effort. 

1.2. The ODE formula 

Let us denote by Q the order of the ODE we are looking for and by D the degree 
of the polynomials in front of the derivatives (we write the ODE in the homogeneous 
derivative x -^). We must then have [Q + 1){D + 1) terms in the series in order to 
determine the unknown polynomial coefficients. If an ODE exists, it appears that the 
number of terms actually necessary for the ODE to be obtained is given by 

N = {Q + l){D + l) -/, (11) 

where / is a positive integer and indicates the number of ODE-solutions to the linear 
system of equations for the polynomial coefficients. 
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From empirical observation, we have seen [7] that N is also given, linearly in 
terms of Q and D, by 

N ^ d-Q +q- D -C. (12) 
While Q and D are the order and the degree, respectively, of any non-minimal order 
ODE that we choose to look for, the parameters d, q and C depend on the series we 
are working with. In all the cases we have considered, we have found that q is the 
order of the minimal order ODE and d is the number of singularities (counted with 
multiplicity) excluding any apparent singularities and the singular point x = 0. The 
parameter C was shown in [8] to be in an exact relationship with the degree Dapp of 
the apparent polynomial of the minimal order ODE 

i^app = id-l){q-l) -C -1. (13) 

Note that there are many ODEs that annihilate a given series. Among all these 
ODEs, there is a unique one of minimal order. In our calculations we have seen that 
it is easier to produce ODEs, which are not of minimal order |11| . in the sense that 
fewer terms are needed to obtain these ODEs compared to what is required to obtain 
the minimal order ODE. Even more importantly for computational purposes, there is 
a non-minimal order ODE that requires the minimum number of terms in order to be 
obtained. 

Next we demonstrate how we use the ODE formula to optimize our calculations, 
i.e. generate just the necessary number of terms in the series. From (jlHllSI) . the 
parameter D is given as: 

D = d-l+ £^ER±A (14) 
Q-q+l ^ ' 

and this must be a positive integer. The parameters / and Q are integers with the 
constraints / > 1 and Q > Q. It is a simple calculation to run through the integers / 
and Q resulting in a positive integer D. For each such triplet (Qoj fa) the number 
■^0 = (Qo + l)(^o + 1) ^ /o is the number of terms in the series required to obtain /o 
ODEs of order Qo and degree Do- Among all these A'o there is a minimum. We call 
the corresponding ODE the "optimal ODE". To obtain the ODE for other primes, it 
is thus only necessary to generate the minimum number of series terms. 
For instance, for x^^\ the ODE formula reads 

N = 12 -Q +33-1) - 887 = (Q + 1){D + 1) - /. (15) 
The optimal ODE, i.e., the ODE that requires the minimum number of terms in the 
series has the triplet {Qo, Dq, Jq) — (56,129,8) which corresponds to the minimum 
number iVo = 7402. Note that the minimal order ODE has the triplet (33, 1456, 1) 
and requires 49537 series terms. 

The minimum number of terms Nq is implicitly given by the ODE formula (jl2[) . 
Plugging the parameter D given in (fTi|) iniV = (Q + l)(-D + l) — /, one obtains 



N= {Q + l)d+ i?app + ^^n-'^l\'^ - (16) 



We can view iV as a continuous function of Q and / and we find that it has two 
extremums when dN/dQ ~ 0. For the positive extremum one has 

Qo = q-l + ^ ^iDapp + f)qd, (17) 



1 



Do = d-l + - ^{Dapp + f)qd, (18) 



Nq^ qd + Dapp + 2 JiDapp + f)qd (19) 
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For the example of x'^' considered above, one obtains (with / = 1) 

Qo ^ 57.20, Do ~ 125.97, Nq ~ 7388.09. (20) 

The gain in the number of terms is already very significant for Q = q + 1 and can 
be measured by the discrete derivative of the hyperbola N{Q) given in (fT6|) . Since we 
should compute over the integers, it is easier to compute the difference of (D+l) (Q+l) 
evaluated at the points Q = q and Q = q + l. At the order Q = q, from one 
obtains D{Q = q) = d — 1 + Dapp + /i, where /i is a positive integer. At the order 
Q = q+l, one has D{Q = q+l)=d~l+ (-Dapp + /2)/2, where /2 is a positive 
integer with same parity as -Dapp- The gain in the number of terms is 

ANiq,q + l)= -d +qh +^{D,pp- f2)q +{fi- f2). (21) 

For x^^\ ^nd with the values fi = 1 and /2 = 2 (since Dapp — 1384 is even), the 
"saving" in the number of terms is 22736 to be compared with the 49537 terms needed 
to obtain the minimal order ODE (i.e. Q = q). As Q increases, one approches the 
minimum of the hyperbola p6)) which is A^o = 7388.09 (with / = 1). Over the 
integers the minimum is 7402 obtained with f — 8. This process can be repeated by 
computing A N{q, q + 2) and in this case D^pp + fz should be multiple of 3. 

As can be seen from the "discrete" derivative (|2T|) . the degree of the apparent 
polynomial is crucial. For ODEs with no apparent singularities the minimal order 
ODE is the optimal ODE. In this case, the hyperbola N{Q) can still have a minimum 
that is not in the integers. 

Note that we may define a minimal degree ODE, i.e. the ODE that has D ~ d 
meaning that there is no singularities other than the "true" singularities of the minimal 
order ODE (no apparent and no spurious singularitie^) . The order of this minimal 
degree ODE is (see p^ ) 

Q = q+ Dapp + / - 1, (22) 

giving for x^^-*, the order Q — 1417 and 103513 as the number of terms (the minimum 
/ being 1). Note that this minimal degree ODE is useless for our computational 
purposes. 

In this paper all of these types of modular calculations and approaches have 
been applied to x'^'- Section 2 shows the computational details (timing, ...) for the 
generation of the first series and the first ODEs, modulo a prime, from which we infer 
the optimal length of the series to be generated for other primes. In Section 3, we 
report on the ODE annihilating x*-^'' and on the ODE annihilating the corresponding 
"depleted" series. The singularities and local exponents confirm the results obtained 
from a diff-Pade analysis and given in a previous paper [7j. In Section 4, the program 
of factorization developed for x'''^'' is used to factorize as far as possible the differential 
operator corresponding to the ODE of x^^''- We will see that our conjecture [8l fTT] 
on the factorization structure of the x^"'* holds for n = 6. Some right factors in 
the differential operator for x*"^"* are obtained in exact arithmetics. Section 5 is the 
conclusion. 

II If we denote by Lq the minimal order differential operator, the non-minimal order differential 
operator Lq^o (with Q > q and D > d) has D — d singularities which are spurious with respect 
to Lq. The spurious singularities are the ones of the operator Lq^^q occurring in the factorization 
Lq,d = Lq-q ■ Lq. 
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2. The series of x'^^ modulo a prime 

As shown in [7] the calculation of a series for x^^^ is a problem with computational 
complexity 0{N'^ In. N). Note that x^^"'' is an even function in w and we therefore 
generally work with a series in the variable x — , though the series for x*-^-* is still 
calculated in the w variable. In Table [T] we have listed a summary of results for the 
formula (11) for various series with new results for x*-®-* added. In [7] we gave a rough 
estimate of the number of terms required to obtain the ODE for x^^'' and thought 
this beyond our computational resources. However, upon closer inspection of Table [T] 
one observes that the minimum number of terms required to find the ODE in x for 
^(2n)^ Yi—12 (or 'I'^"'') is always smaller than the number of terms required for 
^(2n-i) (-Qj. xhis also holds for the combination 6x^"^ ~{n-2) x^"~^^- It is 

reasonable to expect that this would be true for x*-^-* as well. In particular this would 
mean that the number of terms required to find the ODE for Gx^®-* ^ 4x^^'' should be 
smaller than the 6400 or so terms needed to find the ODE for Gx^^-* — Sx^^-*. There 
is of course no way of knowing whether or not this line of reasoning is correct. In 
particular we would have liked to further reduce the number of terms to be calculated 
(one can for instance note that the number of terms required to find the optimal ODE 
for is some 10-20% less than the number of terms required to find the 

optimal ODE for x'^""^'' or ^\ respectively), but since finding the ODE for the 

first time is a hit-or-miss proposition we naturally wanted to ensure, to the greatest 
extent possible, that we had enough terms to find the ODE for 6x'^' — 4x'-^''. For this 
reason it was decided to generate a series to order 6500 in x (13000 in w) for x^^-* with 
the firm hope that this would suffice to find the optimal ODE for at least Gx^^-* — 4x*-'*-' 
(in fact it is also enough terms to find the optimal ODE for x^^-* itself). 

Table 1. Summary of results for various series. The last three columns are the 
data for the optimal ODE. The "l?^' series are the model integrals |6]. 
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N = d- 
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In [7] the calculation of x^^'' to 10000 terms required some 17000 CPU hours on 
an SGI Altrix cluster with l.GGHz Itanium2 processors. Given that the algorithms 
for x^^^ and x*-^^ has the same computational complexity this would indicate that the 
time required to calculate the series for x*-^"* to 13000 terms in w would be at least 
50000 CPU hours (the algorithm for x*"^"* has a slightly larger pre-factor than that 
for x^^^)- III fact it turned out that almost 65000 CPU hours was required and this 
calculation was performed over a six months period. 

The series to order 6500 was calculated modulo the prime 32749. As in jS\ we 
want to factorise various differential operators and reconstruct the right-most factors 
exactly using the results from several primes. We thus need to reduce as much as 
possible the length of the series by identifying some right factors. 

As we detail in the following section the optimal ODE for x^^^ can be obtained 
with less than 6300 terms while the optimal ODE for the combination 6x^^^ — 4x'-^'' 
requires 'just' over 5400 terms. Furthermore we find (using our series modulo a single 
prime) that x^^'* is a solution of this ODE and that one can simplify further by 
considering the linear combination — x'®-* — ^x'"^^ + ^X^^'' whose optimal ODE 
requires a little more than 5100 terms. 

The ODE for x^^'' has as the lower derivative, meaning that Ci + C2 a; is a 
solution (ci and C2 are constants) . Checking that ci + C2 x is still a solution of the 
ODE for ^'^^^ and producing the series ^j^$'^^^(x), we arrive at a series whose minimal 
ODE requires a little less than 5000 terms. 

We therefore calculated a further two series to order 5000 modulo the primes 
32719 and 32717. These calculations required an additional 45000 hours of CPU time. 
Using the factorisation procedure detailed in Section 4, we found a factor of order 3, 
X3, which right divides the differential operator for ^^$'^^■'(2;), and we managed to 
reconstruct X3 in exact arithmetic using 3 primefl. Applying X3, that is form the 

series ^"3 (^^^'I''^''^ 7 results in a series whose optimal ODE requires less than 4800 
terms. At about the same time as these developments took place a new system was 
installed by National Computational Infrastructrure (NCI) whose National Facility 
provides the national peak computing facility for Australian researchers. This new 
system is an SGI XE cluster using quad-core 3.0GHz Intel Harpertown cpus. Our 
code runs about 40% faster (takes about 0.6 times the time) on this facility when 
compared to the Altix cluster and a calculation of a series to order 4800 takes about 
11000 CPU hours per prime. We calculated series to this order for a further 6 primes, 
namely, 32713, 32707, 32693, 32687, 32653 and 32647 (some of these calculations were 
performed on the facilities of the Victorian Partnership for Advanced Computing using 
a cluster with AMD Barcelona 2.3GHz quad core processors). 

3. Fuchsian differential equation for x*^®' 

From the x^^^ series modulo a prime, we obtained various ODEs which have the ODE 
formula 

43Q-I-52D -1121 (Q + 1)(D + 1) -/, (23) 

thus showing that the ODE for x^^'' is of minimal order 52. We denote by ^52 the 
corresponding linear differential operator. 

^ X3 is equivalent to the differential operator L3 given in this paper. 
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The polynomial in front of the highest derivative and carrying the singularities 
of L52 (i.e. the ODE of x'^^) reads 



(1 - 16 xy'' • (1 ~ 25 x) • (1 - 9 a;) • (1 - x) • (1 - 4 x)^ • (1 - 8 x) 

X (1 - x + mx^) ■ (1 - 10a; + 29a;^) -Papp, (24) 

where -Papp is a polynomial whose roots are apparent singularities. Even though we 
have not computed the minimal order ODE, from (|13p . we can infer that the degree 
of Papp is -Dapp = 1020. All the singularities agree with the ones found in 7J from a 
diff-Pade analysis and we have confirmation that (1 — 8x) is the only singularity not 
predicted [5] by the simplified integrals . 

Furthermore, using the exact (modulo a prime) ODE, we can confirm the local 
exponents computed from a diff-Pade analysis in 7] for all singularities except those 
at x = 0,1/16 and a; = 00, which are correct but incomplete. The complete set of 
local exponentf|3 at these latter points read: 

x^o, p= -l,-l/2,o^l/2,l^3/2,2^3^4^5^6^7^8^9^ 

10', 11, •••,17, 

x = l/16, p= -2,-7/4,-3/2,-5/4,-13,-1/2,06,1/2,14,24,33,43, 

53,6^7^8^9^lo, •••,21, 

a; = 00, p = -1^,-1/2, 03, 1/2^ 1^, 3/2^ 2, 5/2^, 3, 7/2^, 4, 9/2^, 

ll/2^ l3/2^ l5/2^ l7/2^ 19/2^, 21/2^, 23/2^, 25/2, 27/2, 

29/2, 31/2, 33/2, 35/2, 37/2, 19. 

Having obtained the ODE formula one can see that the minimal order ODE 
requires 56391 terms (plug Q = g = 52, d = 43, Dapp — 1020 and f — I into HH)). 
And it is a simple calculation, (see paragraph after ((M])) to obtain the number of 
terms necessary for the optimal ODE. This corresponds to Q = 84, D = 73, / = 3 
and N — 6287 terms. If we had to produce the optimal ODE for x'^-* for other primes 
it is 6290 series coefficients that should be generated. 

As mentionned in the previous section, our conjecture that the x^"^ satisfy (with 
an-2 = (n-2)/6) 

X'"^ - ««-2 • X'""'^ + Pn-4 ■ X'""'^ + • • • + <i>("\ (25) 

is also verified. For the series 

$(6) ^ ^(6) _2~(4) ^1^(2)^ (26) 

we obtain non-minimal order ODEs from which we infer the ODE formula 

39 Q -I- 46 £» - 861 = (Q + 1){D + 1) -f, (27) 

showing that the minimal order is 46 with an apparent polynomial (see (|13p ) of 
degree Dapp = 848. The minimal order ODE for <i>(6) requires the generation of 
41736 coefficients series, while the optimal ODE requires 5120 terms corresponding to 
Q = 79, D = 63. It is interesting to see that the required number of terms decreases 
sharply from 41736 (for the minimal order ODE Q — q — 46) to 22272 for the non- 
minimal order ODE Q = q + 1 = 47. The gain AA^(46, 47) = 19464 terms is given 
by ((2T|l for g = 46, d = 39 /i = 1 and /2 = 2 since Dapp — 848 is even. The gain 
AiV(46,48) = 25958. 

t The notation is 0^ for 0, 0, and 9/2^ for 9/2, 9/2, etc. 
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Denoting by L^q the differential operator corresponding to and recalling [TT] 
the differential operator Lio corresponding to^^^\ one sees from (PE)) that the 
differential operator for x^^-* has the "direct sunull decomposition" 

L52 — Lio ® L46. (28) 

The sum of the orders of the differential operators L^q and Lio is larger than 52, 
indicating that a common factor, namely an order four differential operator, occurs 
at the right of both L46 and Lio- The solutions of this order four ODE have been 
given in eqs. (31-33) and eq. (43) of|llj. The differential operator (that we denote by 
L^f^) is given in eq. (42) of [TT] as a product of four order one differential operators. 
Since the expressions for these differential operators were not written in we give, 
for the sake of completeness, in [Appendix A| the full factorization of the differential 
operator L^"*^ . 

Furthermore, we note that in the ODE for x^^"> the derivatives of order zero 
and one are missing (the corresponding differential operator has as the lowest 
derivative^) . The constant and the degree one polynomial x are solutions of L52. The 
constant is a solution of the common factor L^^"^ , but the degree one polynomial x is 
not a solution of Lio and thus should occur in L46. 

We thus have an order five differential operator that right divides L46 

L,= Dl® Lf = (^D^ - ® Lf . (29) 

We now turn to the factorization modulo a prime of the differential operator L46 
keeping in mind that L5 is a right factor. 



4. Factorization modulo a prime of the differential operator L46 

The local exponents at the singularities of the ODE of allows us to easily track 
the factors carrying the various singular behaviours. What we mean is the following. 
The local exponents for the ODE of at for instance w = 0, are all integers. 

Producing the series having the highest exponent, we obtain either the full ODE or a 
right factor. If the series with the highest exponent yields the full ODE then in order 
to obtain a right factor we have to look at the ODEs corresponding to combinations 
of series involving both the highest and the next highest exponent as explained and 
done in [8]. 

For and a.t x — 0, we have two types of local exponents, integer and half- 
integer ones. We thus have a "partition" of the solutions to the full ODE. In other 
words we have "two highest exponents"!^ and it is therefore more likely that we can 
avoid using combination series. 

The ODE for corresponding to L4e has at a; = the local exponents 

p = - 1^-1/2, o^ 1/2, l^ 3/2, 2^3^4^5^6^ 7^8^9^lo^ 11, 

12,13 (30) 

we have then two "highest exponents", p = 13 and p = 3/2. This means that we can 
produce both the series 2;'' (1 -I- • • •) and see whether or not either of these gives rise to 

% Recall that the differential operator for x^^' is a factor in the direct sum of Lio. 
I The notation Dx is 

§ Note that for x'^) , other singularities than w = have half- and fourth-integers exponents. There 
was no need in [8] to use the procedure presented here. 
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a right factor. If so we may not need to resort to the combination method presented 
in Section 4 of [5] . 

At the singularity x — 1/16, the local exponents are 

p= -2, -7/4, -3/2, -5/4, -1^,-1/2, 0^, 1/2, 1^, 23,3^, 4^, 5^, 6^, 

7',8,9,---, 19 

and we have three "highest exponents", p — —5/4, p ~ 1/2 and p ~ 19. 

At the singularity x = cx), there are two "highest exponents" , p — A and p — 33/2 
since the local exponents are 

p = - l^-l/2^0^1/2^, 1^3/2^2, 5/2^3, 7/2^4,9/2^11/2^ 
13/2^ 15/2^, 17/2^, 19/2^, 21/2, ••• , 33/2 

Before we proceed, we introduce the notation L^q — 0„2 • 0„3 , with A6 — ^2 +ri3, 
which we use to indicate that the operator L^q factorizes into two operators of orders 
712 and ria, respectively. Only when a differential operator is definitive do we give it a 
label other than the O. 

Let us begin by the conjecture [8J that L^q has a left-most operator of order six 
which is the symmetric fifth power oi Le- Solutions to the symmetric power oi Le are 
polynomials of homogeneous degrees in the elliptic integrals with the coefficients of the 
combination being rationals. The solutions carrying the half- integer exponents should 
therefore be those of an operator occurring necessarily at the right of L46. So from the 
two "highest exponents" p = 13 and p — 3/2 at a; = 0, we need only obtain the ODE 
of the unique series x^^"^ (1 -I- • • •). Indeed, acting by L^q on the series x^^"^ • (1 -I- • • •) 
produces a series annihilated by an order 40 ODE, leading to the factorization 

LiQ ~ Lq ■ O40. (31) 

When we shift L46 to a; = 1/16 and act on i^^/*^ • (1 -f • • •), with t — x — 1/16, 
we obtain an order five ODE leading to: 

L46 = 041-05- (32) 
Shifting to a; = 1/16 and acting on t^/^ • (1 + • • •) produces: 

^46 = Oae-Oio- (33) 
Shifting Lid to x — oo and acting on i'* • (1 + • • •), with t = 1/a;, gives: 

^46 = 033-Ol3. (34) 

Some factors are common to these three factorizations. Shifting the ODE back to 
a; = and carrying out our factorization procedure [8], one obtains (some final 
labelling is given) 

L46 = 041- O5 = Oii-U-L2, (35) 
Li& = O36 • Oio = O36 ■ Oi • L4 • ^3 ■ -^2, (36) 

i46 = O33 • Ol3 = O33 • Ol • O3 • i4 • i3 • L2. (37) 

The order one differential operator Oi in the last factorization is equivalent to an 
order one differential operator occurring in the Z5 of (|^ . The product O3 • L4 • Z3 • L2 
can be expressed as a direct sum 

O3 • i4 • ^3 • ^2 = ^3 © [U ■ h ■ L2) . (38) 
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Collecting the results given in the factorizations ([31]) and ([57)1 with (|38p. and 
keeping in mind the right factor (j29p . one obtains 

L4Q = Lg • L23 ■ Ln (39) 

with 

Ln = L5 © L3 © (i4 • L3 ■ L2) , (40) 



1 



L5 = (^i?,-^j©Lf. (41) 

Having obtained all these differential operators, a final check is performed by 
acting on by the corresponding ODEs in the order given in ([39]) and doing this 
we do indeed get zero. 

4.I. The differential operator Lq 

The sixth order linear differential operator Lq is the one that we conjectured [8 should 
annihilate a homogeneous polynomial of the complete elliptic integrals E and K of 
(homogeneous) degree five. It should then be irreducible. The local exponents at the 
origin of the linear ODE corresponding to Lq are 

a;= 0, p = -12, -11, -8, -5, -4, (42) 

Plugging a generic series Cn x"" into the linear ODE fixes all the coefficients with 
the exception of the coefficient cq. The "survival" of a single coefficient is a particular 
feature of an irreducible factor with one non-logarithmic solution. The differential 
operator Lq being a symmetric power of Le means that its solution is a polynomial 
in E and K defined as 

K = 2F1 ([1/2, 1/2], [1], 16a;) , E = 2F1 ([1/2,-1/2], [1], 16x) . (43) 

The ODE corresponding to Lq should only have singularities at x = 0,1/16 and 
X = 00, and this is indeed the case. The local exponents at a; = 1/16 are 

X ^ 1/16, p = -48^, -47, -44, -40, 0. (44) 

The local exponents at a: = and x = 1/16 suggest the following ansatz to be 
plugged into the linear ODE (of Lq): 

.i^.(l-16.)^s -E^^--^-W-^^-^- (45) 

The polynomials P^^i_i{x) can be determined numerically and the solution (analytical 
at a; = 0) of the ODE corresponding to Lq is 

+ (1 - IQxf P3,2 -K^E"^ + {1- IQx) P2,3 • K"^ 
+ PiA-KE'^ +Po.5-E'^^ 



The polynomials P5_i,i(x) with coefficients known modulo a prime, are of degree 
respectively, 111, 112, 113, 113, 113 and 113. As conjectured the linear differential 
operator Lq is thus equivalent to the symmetric fifth power of Le- 
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4-2. The differential operator Ln 

The differential operator Ln has in its decomposition the differential operator is 
which is known exactly. The solutions of are the degree one polynomial x and the 
four solutions of l'^'' given in [IT]. As for the other factors of L17, i.e. L2, £3, L3 and 
Li, we have been able to express all of them in exact arithmetics. 

To express a differential operator in exact arithmetics the straightforward 
approach is to rationally reconstruct the differential operator using several modulo 
prime calculations. However, an alternative would be to reconstruct the solutions to 
the differential operator if they are known. This is what we have done for L2 and L3. 

The singularities of the ODEs corresponding to L2 and L3 are only x ~ 0, 
X = 1/16 and x — 00. It is therefore reasonable to assume that the solutions can 
be expressed as polynomials in K{x) and E{x). 

For the ODE corresponding to L2, the solution (analytical at x = 0) written in 
terms of is: 

sol(L2) = (^^-2) X^'^- (46) 

Written in this way, it is easy to recognize the coefficients in exact arithmetics with 
only two primes. The differential operator L2 is thus: 

x-[l~lox) a::-(l — loa;) 

For the third order differential operator L3, we assumed that it is equivalent to 
a symmetric square oi Le- Indeed, the solution (analytical at x = 0) written also in 
terms of (x^^^) , appears as: 



sol(i3) = i • f X (1 - 16xf (16x - 3) • 4^ + (1 - 16x) (64x2 _ 44^; + 9) . A 
X V dx^ dx 

-8(l-8x)(16x + 9)) (x'-^^y ■ (48) 

Here also, two primes are more than sufficient to recognize the coefficients. The 
differential operator L3, in exact arithmetics, reads 

L, = Dl +PlDl +PlD^ (49) 

P3 P3 P3 

with: 

P3 = .x^ • (1 - 16 x)2 (-81 + 1986 X - 17056 .x^ + 34304 x^ + 8192 x"*) , 
P2 = 2. x^ • (1 - 16 x) (2247- 46496 x + 357888x2- 565248 x^- 65536 x"*), 
Pi = 6 (27 - 942 X + 11152 x^ - 101632 x^ + 372736 x* - 65536 x^) , 
Pa = 12 (9 - 308 X- 6208x2 - 101376x3- 49152 x^). 

We have not been able to find the solution of the ODE coresponding to L3. 
The rational reconstruction has been done on the differential operator itself (see 
Appendix B I . Rationally reconstructed, the differential operator L3 reads 

L3 = Dl +^^Dl +'-lD^ (50) 

93 93 93 
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with: 

93 = x2.(l-4x)(l-16a;f Q3, 

03= - 8 + 252 a; - 1678 + 3607 + 4352 x^ 

92 2 a; • (1 - 16 a;)^ (-12 + 1172 a; - 30499 x'^ + 252146 

- 872579 + 770128 x^ + 1 183744 x^) , 

qi = 4 (1 - 16 a;) + 185 x - 28373 x^ + 689440 a;^ - 5128290 x'^ 

+ 16119599 x^ - 13139200 x^ - 17825792 , 
qo = 4 (^-294 + 9469 x + 84480 x^ - 4652220 x^ + 33948640 x^ 

- 97687536 x^ + 89128960 + 74981376 x^) . 

AU the calculations on the previous differential operators have been done with the 
two primes 32749 and 32719. For the differential operator L4 we need more primes. 
The differential operator L4 has the form 

Li = x^ ■ (1-16 x)^{l- Ax) (1-8 x) QiP^^^'' -Di 

+a:2 . (1 - 16 a:)3 QI ■ DI + x [1 - 16 xf Qi • 

+ il-16x)Q;Pt'^ -D, +P^''\ (51) 

where Q3 is the apparent polynomial of L3 in (|50p and pj^^^ are polynomials in a; of 
degree n. To perform the rational reconstruction of the polynomials Pj^^\ we had 
to generate the series for for another seven primes, then obtain the optimal 
ODEs and factorize the differential operators L46 for each prime. After the rational 
reconstruction was completed successfully the resulting differential operator L4 was 
checked against the local exponents and the conditions on the apparent singularities. 
The polynomials P^ "^ are given in exact arithmetics in [Appendix C[ 

Note that we have also checked that these rationally reconstructed differential 
operators are globally nilpotent as they should be. 

4-3. The differential operator L23 

The differential operator L23 has the ODE formula 

21Q +23P» + 1360 = (0 + l)(i:' + l) -/, (52) 
and at X = 0, the local exponents read: 

p = - 25, -24, -23^ -22^, -21^, -20^ -19, -18, -17^, -16, 1, 2, 3, 
4,5,6,-47/2,-45/2 

We can use the same method as before in order to factorize i23- By producing the 
series with the highest local exponents p = 6 and p = —45/2, we obtained the full 
ODE for each series, i.e. an ODE formula compatible with the minimal order 23. 
The singularities of the linear ODE corresponding to L23 are (besides x = 0): 

(1 - 16a:) (1 - 4a;) (1 - a;) (1 - 9a:) (1 - 25a:) (1 - 10a: + 29 a;^) (1 - a; + 16^^) 

We may then shift the ODE corresponding to P23 to a singular point other than a; = 0, 
produce the series of the highest exponent and see whether this gives an ODE of order 



High order Fuchsian ODE 



14 



less than 23. At x — 1/16, the series of the highest exponent p — 11 produced the 
full ODE. Likewise, at other points and exponents such as (x = 1/4, p = —41/2), 
(x = 1/9, p = -47/2), {x = 1/25, p = -63/2), {x = I, p = -47/2) and 
(x = CO, p — —38, —47/2), the series give rise to the full ODE. 

Next we show how the local structure of solutions appear around x — 0. We 
introduce the notation [x^] to mean that the series begins as x^ ■ [const. +•••). The 
results of our computations are the following. Two sets of five solutions can be written 
as (with fc = 1, 2) 

[x*^] In(x)'* + [x-21] ln(x)3 + [x'^^] ln(x)2 + [x'^^] ln(x) + [x'^^], 

[X*^] Mxf + [X~21] 1^(^)2 _^ j^-22] ^j^(^) ^ j^-24]^ 

[x'^] ln(x)2 + [x-21] ln(x) + [x^^*]^ 

[x'^] ln(x) + [x-^i] and [x% (53) 
Three sets of three solutions can be written as (with k — 3, 4, 5) 
[x'^] ln(x)^ + [x-^i] ln(x) + [x"^'^], 

[x*^] ln(x) + [x-21] and [x^^]. (54) 
Two solutions can be written as 

[x^] ln(x) + [x] and [x^] (55) 
Finally there are two non-logarithmic solutions behaving as x"**"^/^ • (1 + • • •) and 

a;-45/2.(i 

Besides the series x'' • (1 -!-•••) with (p = 6 and p — —45/2) that have given the 
full ODE, we may even try the non ambiguous solutions such [x^] in front of ln(x)^ 
and [x^] in front of ln(x)^. But these series produce the full ODE. 

As is the case with the twelfth order differential operator L12 occurring in x^^\ we 
have no final conclusion as to whether or not I/23 is reducible, and without performing 
the factorization based on the combination method presented in Section 4 of [8 we 
do not expect to be able to reach any such conclusion. The representative optimal 
ODE of L23 used in the calculations is of order 67, making the computational time 
obstruction more severe than what we faced with the twelfth order differential operator 
occurring [8] in x^^\ 

4.4- Summary 

Let us now summarize our results. The linear differential operator L46, corresponding 
to «>(^) = x^^^ - fx'"' + ^x'^' can be written as 

L46 — Lq ■ L23 ■ Ln, (56) 



with 



i?. - ^ j © ^3 ® (Li ■ U ■ L2) (57) 



The order seventeen linear differential operator Ln contains only the singularities of 
the linear ODE corresponding to Liq (the operator for x^''') plus the "new"|§| singularity 
X = 1/8. The singularity x = 1/8 occurs only in the fourth order linear differential 

§ It is "new" with respect to what we obtained from the ^^^^ integrals and our Landau singularity 
analysis [7]. 
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operator L4. The third order differential operator L3 is responsible for the p — —5/4, 
p = —7/4 singular behavior around the (anti-)ferromagnetic point x = 1/16. 

Comparing the results of x^^^ with those of x^^\ x^^^ ^-iid x^^^ we note that 
our conjecture still holds: for a given x^""-* there is an order n differential operator 
equivalent to the (n — l)-th symmetric power of at left of the depleted differential 
operators, corresponding to the linear combinations x^^'' ~ B^^^''' X^^"* ~ f^^^'*' 
X^'^ - ix^^^ + TMX(^) and now x^^) - Ix^^) + 

For a given x'"' and once the "contributions" of lower terms (;^("~2fc)^ ^ _ 
1, 2, • • •) have been substracted, the ODE of the "depleted" series still contains some 
factors occurring in the ODE of the lower terms (x'-""^'^')- For x^^^ we have that the 
differential operator Z2 • A^i, which occurs in the ODE of x^^\ continues to be a right 
factor in the ODE of x'^^ — f X^^^ + i^X^^^- For x^^\ we have that the differential 

operator L^^\ which occurs in the ODE of x''^\ continues to be a right factor in the 
ODE ofx('^)-|x(^) + A^(2), 

As was the case for x*-^-* with the differential operators of order two and three {F2 
and i^s), we similarly have for x''^''j the emergence of two differential operators of order 
three and four (L3 and L4), which are globally nilpotent and for which we have no 
solutions. We may imagine that all these ODEs have solutions in terms (of symmetric 
power) of hypergeometric functions (with pull-back) as we succeeded to show |10) for 
Z2. Providing these solutions in terms of modular forms is clearly a challenge. 

Similarly to the twelfth order differential operator L12 occurring in x^^\ we faced 
with the differential operator L23 the same obstruction to its potential factorization, 
namely prohibitive computational times. 

5. Conclusion 

We have calculated, modulo a prime, a long series for the six-particle contribution 
x'^' to the magnetic susceptibility of the square lattice Ising model. This series has 
been used to obtain the Fuschian differential equation that annihilates x'^-* • 

The method of factorization [8] previously used for x^^'' is applied to the 
differential operator L52 of x^^''- With the ODE known modulo a single prime, we 
have been able to go, as far as the computational ressources allow, in the factorization 
of the corresponding differential operator. 

We have found several remarkable results. The factorization structure of L52 
generalizes what we have found for the linear differential operators of x*''^'' ; x'''''' and 
X^^^- In particular, we found in x*'^^ the occurrence of the term x^^'' but also the lower 
term x^^^ leading to the differential operator L46 corresponding to the "depleted" 
series ^'•''^ = x'^^^ - | x^**^ + ^X^^^- The left-most factor Lq of L^q is a sixth order 
operator equivalent to the symmetric fifth power of the second order operator Le 
corresponding to complete elliptic integrals of the first (or second) kind. We expect 
that this happens for all x*'"\ i-^- we conjecture the occurrence in x^"'' of terms 
proportional to meaning a direct sum structure, and the occurrence of a n— th 

order differential operator that left divides the differential operator corresponding to 
the "depleted" series 1^ of x^"^- 

Some right factors of small order appear in the factorization of We have 

used the previously reported "ODE formula" to optimize our calculations. We have 
generated other series of the minimum number of terms, modulo eight other primes, 
and have obtained the corresponding ODEs and the corresponding factorizations. 
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These nine factorizations have been used to perform a rational reconstruction and 
obtain in exact arithmetics the right factors occurring in L^q. 

Our analysis is lacking the factorization of L23 for which, and similarly to L12 
occurring in x*-^-*, we have no conclusion on whether they are reducible. Even if these 
differential operators are known in exact arithmetics, their factorization remains a 
challenge for the methods implemented in various packages of symbolic calculation. 

The massive calculations performed on x^^"* and x^^^ are at the limit of our 
computational ressources and the next step, namely x^^"* and/or x*'^^ seems to be 
really out of reach. A motivation for obtaining these very high order Fuchsian 
operators is to understand hidden mathematical structures from the factors of these 
operators. In this respect, the main results we have obtained on x''^'' are the order 
three and four operators (L3 and L4) that we succeeded to get in exact arithmetics 
and which are waiting for an elliptic curve mathematical interpretation. Providing 
a mathematical interpretation for all these differential operators in terms of modular 
forms is clearly our next challenge. 

The series and differential operators studied in this paper can be found at |12j . 
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Appendix A. The order four differential operator £4^^ 

The order four differential operator L^^^ is a right factor in Lio the differential operator 
for x*"*^- It is a product of an order one differential operator and an order three 
differential operator that can be written as a direct sum: 

Lf = 41(4>4>^x) (A.i) 

In terms of the variable x — they are: 



-(4) 



768x2 

(1 - 16 x) (1 - 24 a; + 96 x'^ 

1 + 384x2 ^ 2048 
2 X ■ (1- 16x)(l- 48x + 128 x^) ' 



and 



= +2^^. (A.4) 
Pi 



with: 

pi = X • (1 - 16 x) (1 - 4 x) (80 X + 7) (-7 + 96 X - 1152 x^ + 10240 x=^), 

Po = 65536000 x^ - 36536320 x^ + 481280 x"* + 254592 x^ - 24800 x^ + 2149 x - 49. 
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Appendix B. Reconstruction in exact arithmetics of the differential 
operator Z3 

The ODE corresponding to is appears as (where the singularities are easily 
recognized): 

+ x-ix- ^) Pi-D, +x-Po (B.l) 
Id 

The polynomials P^,---, Pq arc of degrees, respectively, 4, 6, 7 and 7 in x. We 
have 27 coefficients (not counting the overall one) to reconstruct. For easy labeling, 
these polynomials are denoted as (P3 is the polynomial whose roots are apparent 
singularities) 

3 6 7 7 

P3 = x^ + ^ttkx'', P2 = ^bkx'', Pi = ^CfcX*^, Po = ^dkx''. 

fe=0 k=0 k=0 k=0 

The indicial exponents obtained with both ODEs (with the two primes 32749 and 
32719) are 

a; = 0, p= -2,0,2, 

X = 00, /5 = 1, 2, 5/2, 

X = 1/16, p = -15/4, -13/4, -1, 

a; = 1/4, p= 0,1,7/2, 

Ps{a)=0, p= 0,1,3. 

By demanding that the ODE corresponding to the almost generic Z3 gives the above 
indicial exponents, leads to some conditions on the imknown coefficients ak,bk, Ck and 
dk- The order of the ODE being 3, we obtain for each singularity a maximum of 
three conditions. This is a maximum, because some exponents are by construction 
automatically satisfied. For instance, at a; = 1/4, we obtain only one condition related 
to the exponent p = 7/2. 

At the singularity x = 0, the indicial equation of L3 gives p = as a root 
automatically satisfied and a polynomial in depending on some of the unknown 
coefficients of L3. By requiring p = —2 and p = 2 as roots of this polynomial, we 
obtain 

60 = A ao, CO = ^ ao (B.2) 

With these values assigned, we require that p = 1,2,5/2 be roots of the indicial 
equation at the singularity x = 00. One then gets 

66 = y, C7 = 16, dr = 5. (B.3) 

Similarly, the indicial equations evaluated at the local exponents for the singularities 
X = 1/16 and x = 1/4 give four equations, fixing (e.g.) the coefficients 64, 65, cq and 
dg in terms of other coefficients. 

Next we turn to the apparent singularies. These are given by the roots of P3. 
Calling a a root of P3 (with unknown 0^), the indicial equation appears with p = 
and p = I as automatically satisfied roots. Requiring p = 3 as root of the indicial 
equation, gives a polynomial in a of degree three. Zeroing each term gives 22 solutions. 
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Discarding all the solutions where a coefRcient from is zero, one is left with five 
solutions. From these solutions, there is only one solution which is acceptable, because 
it matches with the actual values of the coefficients known in prime. This fixes three 
coefficients in terms of the others. 

At this point, we have fixed 12 coefficients among the 27 using only the knowledge 
about the local exponents. The condition on the local exponents at the apparent 
singularities is only necessary, the sufficient condition is the absence of logarithmic 
solutions around the singularity x = a. 

The conditions on the non-occurrence of logarithmic solutions at the apparent 
singularities can be imposed either by requiring the conditions of eq. (A. 8) in [7] to 
be fulfilled or equivalently by zeroing the coefficients in front of the log's in the formal 
solutions of L3 at a. With a generic apparent polynomial, the calculations can be 
cumbersome. So let us fix some coefficients. 

One finds that the ratio — 2ai/ao appears with both primes 32749 and 32719 as 
the number 63. Also for both primes one obtains 4a2/ao = 839, —Sa^/ao — 3607 
and 214 ^^/^^ ^ 

Furthermore, one may compute the (analytical at x = 0) series at both primes in 
the hope that some coefficients will be "simple" enough to be recognized. The series 
with the prime 32749 gives 

+ A8x^ + 1527 X* + 7541 x^ + 3199 x^ + ■ ■ ■ (B.4) 

while with the prime 32719, it reads 

x^ + A8x^ + 1527 x^ + 7571 x^ + 4069 x^ + ■ ■ ■ (B.5) 

We note that the same values occur at orders 3 and 4. These numbers are therefore 
likely to be exact. Also the difference between the values at order 5 is a multiple of the 
difference 32749 — 32719, and similarly at order 6. It is easy to "guess" these values 
as respectively, 48, 1527, 40290 and 952920. Comparing with the series solution of L3 
fixes four coefficients. 

We have then twelve coefficients fixed exactly and nine coefficients fixed by 
reconstruction. The formal solutions of L3 at the apparent singularity a give two 
logarithmic solutions, with leading term, each 

C a'' {x-af lii{x-a), /c = 0, • • • , 3 (B.6) 

where C depends on the remaining non fixed coefficients of L3. We have then eight 
(non-linear) equations for six unknowns to solve. This can be done by rational 
reconstruction and check. 

Appendix C. The differential operator L4 in exact arithmetics 

The degree n polynomials pj"'\x) occurring in the differential operator L4 read: 

P]^*^^ = 28000 - 7854000 x + 873083400 x^ - 54037012120 x^ + 2099285510560 x^ 
-52582954690298 x^ + 766418384173454 x^ - 1305110830870633 
-251264549473230968 x^ + 7727889974481947660 x^ 
-148605250583921845896x1° + 2252938824290334087840 a;" 
-29645475671183771992224 a;i2 -I- 354446803792968575565792 a;" 
-3850023960384577768909952x1" + 36761552740911534545901568 a;i^ 
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-296338746597146803591135232 x^^ + 1953967934450852091348254720 a;^^ 
- 10332892566359614848157876224 a;^^ + 43345424617004971574289235968 x^^ 
-142807225508285034141616963584 x'^° + 359505820412663945726355570688 x'^^ 
-636026962079787427490890252288 2:22 + 616797192523902897669611192320 a;^^ 
+45081769872830521912080728064 x^^ - 724445324775545659452335063040 x^^ 
+521686412421099571093753036800 a;26, 

P^^^^ = -4480000 + 1569568000 x - 238072889600 x^ + 21281848471520 x^ 
-1268595101537120 x'^ + 53555230610961720 x'" - 1640958998875092768 x^ 
+36032181180727162732 a;^ - 511428562675996247108 
+1919885419260765103140 x^ + 129005457127386313373184 x^° 
-4541113747259527374959592 a;" + 96035689755227434986877112 x^'^ 
-1580421468708164786235613784 .t" + 22087897691588508601005658336 x^^ 
-274269909442085751262554453856 x^^ + 3087338965228238905750107987648 x^^ 
-31474922613166692487806647824256 a:^'' 
+286292076483602608978320943481344 a;^^ 
-2277952733740370146287983798312960 a ' 
+15571420858521621122719931928608768 
-90147310596750652057735905075527680 
+437037767842717994841190340774330368 
-1755686044559298411692425577783885824 x' 
+5764646607249819312839063743970148352 a;2^ 
-15099644256008129321411266837095645184 " 
+29988658044590195137583404663925899264 
-39745090934862435542362760545321353216 
+19398167217699147074111209484113149952 
+40447257076217292533523320942836580352 
-84060042791775646091063152898350252032 
+40724840987587942318458896159738953728 
+34088801304111660197683822288919592960 
-34873025538917765121024203000119296000 

P^^^^ = -4480000 + 1569568000 x - 238072889600 a;^ + 21281848471520 x^ 
-1268595101537120 a;* + 53555230610961720 a;^ - 1640958998875092768 a;^ 
+36032181180727162732 .t^ - 511428562675996247108 x^ 
+19198854192607651031402;'' + 1290054571273863133731842:1° 
-4541113747259527374959592 2:" + 96035689755227434986877112x^2 
-1580421468708164786235613784x1^ + 22087897691588508601005658336 x" 
-274269909442085751262554453856 x^^ + 3087338965228238905750107987648 x 
-31474922613166692487806647824256x1^ 



.25 

1x26 

ix27 
ix28 
X29 
lX^° 

iX^' 
)X32 
„33 



High order Fuchsian ODE 20 

+286292076483602608978320943481344 a;^* 
-2277952733740370146287983798312960 a;^'' 
+15571420858521621122719931928608768 
-90147310596750652057735905075527680x21 
+437037767842717994841190340774330368 a;22 
- 1755686044559298411692425577783885824 x"^^ 
+5764646607249819312839063743970148352 a;2^ 
- 15099644256008129321411266837095645184 a;^^ 
+29988658044590195137583404663925899264x26 
-39745090934862435542362760545321353216x2^ 
+19398167217699147074111209484113149952x2^ 
+40447257076217292533523320942836580352x29 
-84060042791775646091063152898350252032 

+40724840987587942318458896159738953728x^1 

+34088801304111660197683822288919592960x^2 

-34873025538917765121024203000119296000 x3^ 

P^^^^ = 202496000- 84671104000 x+ 15961404659200x2 - 1817819283938560x3 
+141042261097575040 x" - 7945786419559994432 x^ + 336970482890735391136 x^ 
-10948102706558839518064 x^ + 272101251799491505044720 x* 
-4990110947182458236154960x9 + 57747165172968723279034760x1° 
+4251375690841730042108460 x" - 19664102412813111595220034000 x^^ 
+586539601535060491103255831780x13- 11648280832868820874871506994648 x" 
+185168754164459407231412940918408x1^ 
-2524027149739792644483439537740736x1*^ 
+30612264160202676427790224166656736x1^ 
-336756368430758251349398256374549440x1^ 
+3374928905004383352843307682288939648x19 
-30594829577694461795851875047251759104x2° 
+247550135999641906395555078053550042624x21 
- 1761694860791556801623390940580862476288 x22 
+10880585300165439414813579169355207311360x23 
-57651463831886251900194559835893548711936x2^ 
+259270510361927197193957311877476593434624x2^ 
-977978052427489585499761822245900827754496x26 
+3043305515555663644471442318841392857612288x2^ 
-7593091629989468917294503828609326603304960x28 
+14317720902933442365662690637880059263713280x29 
-17337418172194871339769688830180251691646976x30 
+3546879809404692840748046019057281136590848x31 



High order Fuchsian ODE 21 

+32904223733304447370184725984806679848419328 a;^^ 

-61070255095717193234874579385327453575577600 a;^^ 

+27287160011587832026533318214423160423448576 x^"* 

+47538188516382446352727572349627507901726720x3^ 

-56445574686008125172119780480189438504206336x3'^ 

-4318904703692797702055795738669690860339200x3^ 

+23615008551819589708322104774634383815475200x3^, 

P^('*3) = -2508800000+ 1313872896000 X- 313495056179200x2 

+45402581315051520x3 - 4502030899704432640 x* + 326696241278915100672 x^ 

-18076764858722283537408 x*^ + 782686127310817603163904 x^ 

-26913654199485748976447296 x^ + 737895426074343351817982240 x^ 

-15941906513987915790530627104 x^" + 258773331815879690900773968400 x" 

-2607962306360230233373492782176 x^^ _ 5952736815704779243433578988240 x" 

+986994067078072761785220512495568 x" 

-27146844884553280870530528088810192 x^^ 

+517853080131584647940304813906843912x^*^ 

-8012481021063135055260920360860291792 x^^ 

+106770798207835443855237151398845884336x1^ 

-1266450623115899739824560105189293118336x1^ 

+13633188914202542686825030295715897195712x2° 

-134327854309114583892390549684219213327360x^1 

+1210064015789594600623288132568732268617728x^2 

-9885950925613173310943030286291198265745664x^3 

+72409258384998425181940033322470418064579584x24 

-469743224562744760675515167582702668512534528x25 

+2668377435339727605940145954082991900173762560x26 

-13130369247854115699188867418934857934469332992x2^ 

+55351059606988527054614355506855140926785847296 x2S 

-197223974876465329508996857363261585516738379776x29 

+582489670248346892198679343375535443002292961280x3" 

-1378041300571967278991550115889940047326999478272x31 

+2424920581794299143009574014342980105306252509184x32 

-2497995923785565959357957923036374193256014020608x33 

-915221239768968177447587513000776938544384966656x34 

+8834509277491743877951520843172281230968074797056x35 

-14739433233907061551935598195219565259254451404800 x3'5 

+5672247850181350880926829983674652540556997033984x3^ 

+16865676565374917674763783610716646680263347142656x38 

-23057876384647717319687179507181219615528382889984x3'' 



High order Fuchsian ODE 

-1650497416603706423024253737913692788809736912896 a;^° 

+17950328578610802277327697770535083294270546771968 a;''^ 
-3724970001243182786619005376755715492471782768640x^2 

-5981413341400069058756778898532874294556360704000 a:^^ 

Pq(^^)/16 = 58841859686400 a;- 123282432000- 13099552866570240 a;^ 
+1817269523720161280 - 176880012691621796864 .t"* 
+12880441893460632329216 .T'"^ - 729910851392566766105088 
+33014141392329879832166784 a;^ - 1210942171584302533599014752 a;* 
+36308844474092544015578885632 - 889129088058672919373638221264 a 
+17508381271590013090109310169040 a;^^ 
-263494886656617518206756373588932 a;^2 



20 



+2493591715504400008185935972185648 a;i3 

+5772652777357046820837948335210000 a;i4 
-880112646375062548999453842320020740 a;^^ 
+23493494316713860255651067234081149257 a;^^ 
-438058417861614862884693737035489286345 a;^^ 
+6643943767863126335261566491851505292189 a;^^ 
-86756699901268061114746560625886582904365 a;^^ 
+1006186625234680761751520210312145980149549 x 
- 10573522222931 154420271493264607253396390520 x^^ 
+101894357518227690884588911318694326634020120 a;22 
-904489874897014837177389619617321458811380360 a;^^ 
+7377025197157259422622822297421365236335307120 a;^^ 
-54857750379533672661182684179932897350723993600x2'^ 
+368077157510764846299472690339090755869412496960 .T^e 
-2203473576836831766402446311571835988588370992640x2'' 
+11638948368194240385082022186232592838207372154880x2^ 
-53640316561668843524196008022033100643589470191616x2^^ 
+213035315257870225008043406258186703365521254907904x^° 
-717584853510007605413068883853037631249102250442752x31 
+2000779126900084461641442809746125018650394950107136x^2 

-4414463297786097513664235192893813927566161255333888x^3 
+6904891435787610921130882736916279097844736823656448x3^ 
-4551724050684467601081502988404586537388373763424256x35 
- 1 1571837065995769727688883612933393577503693010370560 x^^ 

+43604071314966497936511910817544815142484611319201792 x3^ 
--64046695475293378492360343847354456353947960484036608x38 
+17229520899952062417015850756255391466062797246300160x39 
+98913305027317465024954824787190137389180923821424640 x*° 



High order Fuchsian ODE 23 

-145693357979556257119098588624331861246512084740472832 a;-*! 
-1743763200037842518500493452602647084799741447372800 
+159968299464829816606333313819738117801053481636724736 a:'*^ 
-68453133710464189864730237770717937227743579749744640 x^'^ 
-84974525390992986946108353023934616304288806232653824 x^^ 
+41407097440632033071894561954752886956699467613470720 
+28698609854675644415679733396189051258415886630912000 a;^^. 
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